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ABSTRACT 

An analysis of the interaction between a spherical relativistic blast-wave shell and a stationary cloud 
with a spherical cap geometry is performed assuming that the cloud width A c / <C x, where x is the 
distance of the cloud from the GRB explosion center. The interaction is divided into three phases: 
(1) a collision phase with both forward and reverse shocks; (2) a penetration phase when either the 
reverse shock has crossed the shell while the forward shock continues to cross the cloud, or vice versa; 
and (3) an expansion phase when, both shocks having crossed the cloud and shell, the shocked fluid 
expands. Temporally evolving spectral energy distributions (SEDs) are calculated for the problem of 
the interaction of a blast-wave shell with clouds that subtend large and small angles compared with 
the Doppler(-cone) angle Qo = 1/Tq, where Tq is the coasting Lorentz factor. The Lorentz factor 
evolution of the shell/cloud collision is treated in the adiabatic limit. Behavior of the light curves 
and SEDs on, e.g., To, shell- width parameter n, where Ao + ^x/Tq is the blast- wave shell width, 
and properties and locations of the cloud is examined. Short timescale variability (STV) in GRB 
light curves, including ~ 100 keV 7-ray pulses observed with BATSE and delayed ~ 1 keV X-ray 
flares found with Swift, can be explained by emissions from an external shock formed by the GRB 
blast wave colliding with small density inhomogeneities in the "frozen pulse" approximation (rj — ► 0), 
and perhaps in the thin-shell approximation (rj w I /To), but not when 77 m 1. If the frozen-pulse 
approximation is valid, then external shock processes could make the dominant prompt and afterglow 
emissions in GRB light curves, consistent with short delay two-step collapse models for GRBs. 
Subject headings: gamma ray bursts — radiation processes: nonthermal — hydrodynamics — relativity 



1. INTRODUCTION 

Important knowledge about the nature of GRBs comes 
from analyses of burst light curves. BATSE, triggering 
on peak flux over 64, 256, and 1024 ms timescales in the 
50 - 300 keV band, provides the largest database at hard 
X-ray/soft 7-ray energies, amounting to over 2292 GRB s 
reported nearly 8 years into the mission (|Fishmanl ll999). 
and ~ 2700 BATSE GRBs in total. Examining BATSE 
GRB light curves gives the impression that the 7-ray ac- 
tivity, so intense during the first ~ 1 - 100 s, begins to 
weaken with time, all the while displaying pulses with rel- 
atively constant widths. If this evidence requires for its 
explanation a very powerful structured relativistic wind 
from an active central engine that becomes less vigorous 
with time, then the implied GRB explosion mechanism 
is very different than if this data is taken as evidence 
for a single explosive event where the GRB pulses are 
made by external shocks from a weakening blast wave 
that interacts with material in the circumburst medium. 

The GRB afterglow revolution, initiated by Beppo- 
SAX and continued with HETE-2 and INTEGRAL, gave 
for the first time GRB source redshifts and distances and 
therefore apparent isotropic powers and energies. With 
the opening angle of the relativistic jet inferred from 
the achromatic beaming breaks in optical light curves, 
the absolute energy of a GRB is fixed to the uncer- 
tainty in the efficiency of radiation production. The 
tools of standard relativistic blast-wave physics for the 
afterglow external shock emissions allow one to test for 
uniform or wind-formed circumburst media, and deter- 
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mine the e e and eg paramete rs for assumed jet structures 
(|Panaitescu &; Kumar 2002) . Fitting the statistical data 
on the GRB redshift and opening angle distributions can 
determine the intrinsic jet opening angle and whether a 
top- hat model for GRB jet structure gives good agree- 
ment with the data (|Guetta fc P iran 2006]; iLe &: Dermer! 
l2007lh 

A rich new data base is recently opened by Swift 
(|Gehrels et al.ll2004[) . After triggering with BAT in the 
15 - 150 keV range, Swift autonomously slews within 
« 100 s to a burst so that the XRT, with excellent sen- 
sitivity, can get a clear picture of the evolution of the 
0.3 - 10 keV X-ray emission. In large numbers of GRBs, 
there are steep declines in the X-ray emission at the end 
of the prompt X-ray luminous phase hundreds of sec - 
onds after the start of the GRB dTagliaferri et al.ll2005h . 
an extended plateau phase dBarthelmv et al.ll2005h . and 
X-ray flares ([Burrows eLaDl2005), which in some cases 
have durations that are a small fraction of the time since 
the start of the GRB. Extrapolation of the BAT into the 
XRT band gives an almost continuou s X-ray light since 
the GRB trigger (|Q'Brien et al.ll2006h . 

The interpretation of this data is complicated not only 
by the rich detail revealed by Swift but, for statistical 
studies, by its unusual triggering method, which makes 
the a ssignment of a threshold flux problematical ([Band! 
200(| . Because the Swift BAT uses a coded mask from 
which images of the GRB photons can be formed, a lower 
peak flux than pre-Swift telescopes can be identifiable as 
a GRB by using both rate and image triggers. A broader 
dynamic range, from 0.004 s to 26 s, is used for the 
rate trigger. Differences in redshift distributions of Swift 
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GRBs compared with the GRBs detected with BATSE, 
Beppo-SAX and HETE-2, which had comparable trig- 
gering criteria can, however, mainly be understood by 
differen t peak flux sensitivit ies of Swift and pre-Swift de- 
tectors ULe fc Dermerll2007h . 

The X-ray flares in the late prompt and early after- 
glow (t ~ 10 2 s - 10 4 s) observed with the Swift satel- 
lite have been interpreted as evidence for internal shocks 
and refreshed energization of the relati vi stic outflow 
Zhang et all 120061: IChincarini et aI1l2007f ). iLiang et all 



20061 ) and lYamazaki et alj l 20061 ) model the steeply de 



dining phase as a superposition of background external 
shock emission and a curvature pulse produced by shells 
ejected long after the start of the GRB. This model for 
energy dissipation in GRBs and production of the GRB 
light curves is commensurate with the collapsar/failed 
supernova scenario for long-duration GRBs, 1 where the 
evolved iron core of a massive star collapses to a black 
hole. In the collapsar picture, energy dissipation in an 
accretion torus or through the Blandford-Znajek process 
drives a narrow collimated jet through the stellar enve- 
lope. During the early prompt phase, continued activity 
drives successive waves of collimated relativistic plasma 
through the SN shell. Collisions between the shells make 
GRB pulses. Later activity and peculiarities in the light 
curves are attributed to refreshed outflows or internal 
shocks from late outbursts of the engine. Interpreting 
rapid variations in the GRB X-ray light curve as activ- 
ity of the central engine, then the engines in some GRB 
sources must be on f or rest-frame times < 10 4 s after 
the start of the event (|Burrows et alj|2005t iFalcone et al.l 
l2006UBurrows et al.ll2005af) . 

Early calculations of the spectral and temporal behav- 
ior of th e GRB pulses in an internal s hock model were 
made by I Dai gne &: Mochkovitchl ( 1998). Kobavashi et al.l 
(|1997| ). and iGranot et all (119991 1 The efficiency dif- 
ficulties for the internal shock mode l (jKumarl 119991 : 
iBeloborodovl 120001: IMimica et all l2007f) are reduced by 
large Lorentz-factor con trasts between different shells 
(Kobavashi & Sari 2001), which constrains the collision 
radius and vF v pea k energies under th e requirement of 
Thomson thinness (|Guetta et al.|[200ll ). For the inter- 
nal shock/colliding shell model, as well as for the exter- 
nal shock model, efficiency concerns are ameliorated by 
the large hadronic energy content in the prompt phase 
that escapes t o form the ultra-high energy cosmic rays 
(jDerme 3[2003). 

A complete description of the collision between two rel- 
ativistic shells that follows both the hydrodynamics and 
radiation physics needed to calculate emission spectra 
and light curves from different pulses in GRBs, necessary 
to test the internal shell model, h as only been r e cently 
treated, in a series of papers by IMimica et alj (|2004 
2005, 2007). Here we present a complete and systematic 
analysis of the collision between a cold GRB blast wave 
shell and a cloud with a spherical cap geometry in an 
external shock model. We use these results to check the 
widespread claims by the Swift team that central engine 
activity is r esponsible for the erratic X-ray flares from 
GRBs (e.g., lO'Brien et al.|[2006at IFalcone eTail l2006al : 

1 Long-duration GRBs are referred to henceforth as GRBs; here 
we do not consider the short hard GRB class or the low-lumino sity 
GRBs, also argued to form a separate class (Liang et al. 2007). 



iRomano et ail 120061 : [Burrows et al.l I2007D . We dispute 
these claims and argue that an external shock model can 
make the prompt 7-ray pulses and afterglow X-ray flares 
under the condition of no spreading of the cold blast- wave 
fluid shell. If this assumption is allowed, then interest- 
ing implications for the nature of the collapse process 
in GRB sources follow. A short delay two-step collapse 
process, in fact, provides a competing and compelling 
alternative to the collapsar/failed supernova model, as 
discussed below (Section 5). 

The emissivities calculated from the interaction of a 
blast-wave shell and a stationary cloud are integrated 
over three spatial dimensions to obtain spectral fluxes as 
a function of observer time. Only the synchrotron com- 
ponent is considered here (for a recent treatment of the 
syn chrotron-self Comp ton component in external shocks, 
see lGalli fc Piroll2007t ). and other assumptions are made 
to simplify the analysis, for example, a randomly ordered 
magnetic field and isotropic electron distribution in the 
proper frame of the shocked fluid. 

The tec hniques of the standard blast - wave model are 
used (e.g. iSari et all 119981 : iPiranl I1999L 120051: iMeszarosI 
2006), and the results are applied to observations of 
prompt ^100 keV emission pulses and ~ 1 keV X-ray 
flares in GRB light curves. The problem is analyzed 
and broken into three distinct interaction phases, as de- 
scribed in Section 2. Calculations of light curves and 
SEDs from the emissions of external shocks formed by 
the interaction of a GRB blast wave with a stationary 
cloud are presented in Section 3. Parameter are chosen 
to produce short timescale variability (STV; At v /i <C 1, 
where At v is the variability time scale and t is the time 
since the start of the GRB) in the 7-ray lightcurve, which 
is only possible if the blast wave shell width can be ap- 
proximated by a a thin shell (A(a;) < x/Tq) or frozen 
(unspreading) shell. This result is demonstrated analyt- 
ically in Section 4, and illustrated by Monte Carlo simu- 
lations of GRB light curves. We conclude that the X-ray 
flares observed with Swift are compatible with an exter- 
nal shock model if the GRB blast wave shell interacts 
with a clumpy surrounding external medium and does 
not spread transversely. This assumption is discussed in 
Section 5, and is argued to be correct. Implications for 
the nature of the collapse mechanism and jet formation 
in long-duration GRBs are also discussed in Section 5. 
The study is summarized in Section 6. 

The reader who is not interested in the detailed analy- 
sis may avoid Section 2, and refer directly to the numer- 
ical results in Section 3. The principal conclusions of the 
analysis are found in the simplified analytic descriptions 
in Section 4. 

2. ANALYSIS OF THE BLAST- WAVE/CLOUD 
INTERACTION 

Consider a GRB that takes place at redshift z and re- 
leases energy over a characteristic timescale Ao/c repre- 
senting the period of activity of the GRB central engine. 
The corresponding spatial dimension Ao should probably 
be greater than the Schwarzschild radius of the several 
Solar mass black hole formed in the GRB event; thus 
Ao 10 6 cm. Depending on the nature of the central 
engine of GRBs, Ao could range from a fraction of a 
second to hundreds of seconds or longer if, as in the col- 
lapsar model for the GRB prompt phase, the duration 
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Fig. 1. — Illustration (not drawn to scale) of the geometry of 
a spherical, relativistic blast-wave shell interacting with a cloud 
having a spherical cap geometry. The axis through the center of 
the cloud is offset from the observer's line of sight by the angle 
The cloud subtends an angle 6 C ; as measured from the explosion 
center. Emission originates from radiating plasma located at a 
distance x from the center of the GRB explosion, and at angles 8 
and <j> measured with respect to the line-of-sight direction to the 
observer. To approximate a quasi-spherical cloud, the transverse 
and lateral extents of the cloud are set equal. 



of the highly variable X-ray and 7-ray flux in a GRB is 
assumed to reflect the period of activity of the engine. 
Here we consider a GRB engine where the progenitor 
neutron star collapses impulsively to a black hole, so that 
A ~ 10 7 A 7 cm, with A 7 ~ f . 

The apparent isotropic equivalent 7-ray energy re- 
leased by a GRB explosion is written as Eq. The appar- 
ent 20 keV - 2 MeV isotropic rest-frame energy measured 
from GRBs typically range in values from « f 51 - 10 54 
ergs, with a handful of anomalously low energy GRBs 
with energies as low as 10 48 ergs (e.g., GRB 980425; 
iFriedman fe Bloomll2005t iGhirlanda et alj|2007l ). A sig- 
nificant number of GRBs have Eg > 10 54 ergs, so that 
the total GRB energy must be larger, by a factor of at 
least several (and possibly much larger if the radiation 
efficiency is low or ener getic hadrons are formed and es- 
cape, e.g., as neutrons (jAtovan fc Dermerll2003l : iDermerl 
2007)). In this study, jet effects are considered as a re- 
striction on the interaction angle, the jet structure is as- 
sumed to be 'top-hat,' and there is no lateral and very 
little transverse spreading. 

The explosion is assumed to form a fireball with initial 
Lorentz factor (entropy per baryon) denoted by 



r n = 



-En 



»1, 



(1) 



where Mq is the amount of baryonic matter mixed into 
the initial explosion. In this anal ysis, we do not con- 
sider effects of neutron de c oupling dDerishev et al1ll999t 
iBahcall fc Meszarosl [20001: lBeloborodovll2003D . and fur- 
thermore assume that the blast- wave energy in the coast- 
ing stage is carri ed primarily by particle rathe r than field 
energy (see, e.g.. lLvutikov fc Blackman1l200lL for the lat- 
ter case). 

The blast-wave plasma shell reaches its coasting 
Lorentz factor Tq at distance x > To An from the ex- 
plosion. At x > x spr = TqAq, where x spr is the spread- 



ing radius, internal motions within the blast-wave shell 2 
are thought to cause i t to spread s o that A(x) « r/x/Tl 
(jMeszaros et all 119931; lPiranlll999h . where rj < 1. We 
consider the case that the shell experiences little spread- 
ing, so that 77 1 , which has as its asymptotic limit the 
frozen-pulse approximation (77 — > 0). More discussion 
about rj and shell spreading is given in Section 4.4. 

In the stationary frame of the explosion, the shell width 
A(x) is therefore given by 



A(x) S A + 77 



(2) 



The proper number density of the relativistic shell is 
given by 

2.1. Geometry of the Blast Wave/Cloud Interaction 

The interaction event is sketched in Fig. 1. The blast- 
wave shell of width A = A(x ) collides with a uniform 
cloud with density n c /. The cloud is assumed to have 
a spherical cap geometry and to subtend a solid angle 
AQ, c i = 7r# 2 j as measured from the explosion center, and 
it therefore presents a projected area A c i — tt6%Xq to 
the GRB blast wave, where xo is the distance from the 
origin to the inner edge of the cloud. The angle between 
the axis through the cloud center and the line-of-sight to 
the observer is denoted 6*^. The coordinates x = (x, 6, <j>) 
defining the location of the radiating material in the sta- 
tionary (explosion) frame are measured with respect to 
the line-of-sight to the observer, with the angle <fi repre- 
senting the projection of the azimuth on the plane normal 
to the direction to the observer; thus <f> = is defined 
with respect to the projection of the axis through the 
cloud center on this plane (see Fig. la). For simplicity, 
both the shell and cloud are assumed to be composed of 
electron-proton plasma. 

The vF v flux, denoted by f e (t), measured at observer 
time t and at dimensionless photon energy e = hv /m e c 2 , 
is given by 



/«(«) 



2tt 

, , , d(f> I dn I dx x' 2 e'j'(e',CL';x,t') 6'& , 



(4) 

where primes refer to comoving quantities, e' = e z /Si) = 
(1 + z)e/5-D, and the integration is over volume in 
the stationary (explosion) frame jlGranot et al.l 119991 : 
iKumar fc Panaitescd[200a iDerm cr 2004). The Doppler 
factor of the radiating fluid is defined by the expression 



r(i - 



(5) 



where /i = cos 8, T = 



Lorentz factor, and (3c = c\f\ 
differential emissivity 

j'(e / ,n / ;x,t') = 



r(x,t*) is the emitting fluid's 
~ 2 is its speed. 



The 



d£'(n') 



dV'dt'dQ'de' 

2 The relativistic blast-wave shell from the GRB will always be 
referred to as the shell, and the stationary material that is swept 
up by the shell will be referred to as a cloud, even when the cloud 
represents supernova remnant (shell) material. 
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is defined such that j'(e', Q'; x, t')dV dt' d£l' de' is the dif- 
ferential comoving energy d£' of photons with dimcn- 
sionless energies between e' and e' + dt' that are radiated 
during time dt' from comoving volume dV' within the 
solid angle element dil' in the direction ft' defined in the 
comoving fluid frame. Although the directional proper- 
ties of the radiation must be considered in the general 
case of nonthermal synchrotron emission from relativis- 
tic electrons in an ordered magnetic field, or for external 
Compton processes, here we assume that the radiation 
is emitted isotropically in the comoving frame, and as- 
sume a nonthermal synchrotron origin for the keV - MeV 
emission from GRBs. Implicit in this approximation is 
that the magnetic field is randomly ordered and that the 
electrons have an isotropic pitch angle distribution in the 
fluid frame. 

Thus j'(e',Q,';x,t') = j'(e';x,t')/4ir. Because the 
azimuthal dependence no longer appears in either 
j'(e';x,t') or <5d, it is trivial to integrate eq. ((4]) over 
4> under the assumption that the emissivity j'(e'; x,t') — 
j'(e'; x,t'), that is, the emissivity depends only on x and 
t but not on angle (except as determined by the physical 
extent of cloud). This assumption allows us to treat the 
hydrodynamics of the blast- wave/cloud interaction in a 
planar geometry (Section 2.2). We obtain 



fS) = (27T4)" 1 / M\sme\<t>+ 



x dx x 2 e'j'(e';x,t') 5^ 
Jo 



(6) 



noting that the Doppler factor is, in general, dependent 
on location. The term (j> + = cf> + (9,9 c i,9i) represents the 
maximal azimuthal angle subtended by the cloud for re- 
ceived emission from radiating plasma located at an angle 
6 with respect to the observer's line of sight; see Fig. 1. 
From the addition formula for angles in spherical geome- 
try, we have that cos Q c i = cos Q cos Qi + sin 9 sin Qi cos <f>+ . 
Inspection of Fig. 1 reveals that when Qi < 9 c i, that is, 
when the line of sight to the observer intersects the cloud 
volume, then </> + = n when 9 < Q c i — Qi. Thus we have 
two cases, as follows: 



1. 



< Qd 

7T , 



when 



<9 <9, 



b + = { arccos 8 , when Q c \ — 9i < 9 < 9i + 9 c i 

, otherwise; and 

2. Qi > Q c i 

( arccos , when Qi — Q c i < 9 < 9 c i + 9i , 
<f>+ = \ 



(7) 



0, 



where 



9 



otherwise, 



cos u c i — cos a cos ( 



sin u sin ( 



(8) 



The stationary frame time is related to the observer 
time t through 



_ t x cos Q 

tz = TZ T — 



(1 + *) 



(9) 
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Fig. 2. — Phases of the Interaction Event. Top figure illustrates 
the planar geometry and widths, in the stationary frame, of the 
shell and cloud. Panel (a) illustrates the Collision Phase. Phase 
1, with both forward shock (FS) and reverse shock (RS). Panel 
(b) illustrates the Penetration Phase, Phase 2, when either the RS 
has traversed the shell and the FS continues to pass through the 
cloud (left), or when the FS has traversed the cloud and the RS 
continues to pass through the shell (right). Panel (c) illustrates 
the Expansion Phase, Phase 3, where the shock fluid expands and 
the nonthermal particles cool, for the former (left) case of Phase 2. 

where the zero of time in the two frames is defined by 
the start of the GRB explosion. The differential distance 
dx traveled by a relativistic blast wave during the differ- 
ential time elements dt*, dt' and dt x in the stationary, 
comoving, and observer frames, respectively, is given by 

(10) 

Eqs. © and (fTU)) are needed to determine the emissivity 
at assigned values of x and t z . 

2.2. Phases of the Interaction Event 

The complete interaction is treated in a planar geome- 
try, as illustrated in Figs. 1 and 2. The cloud is assumed 
to have uniform density between radii xq and X2 , so that 
the cloud width A c / = X2 — xq. A major simplifying 
assumption of the analysis is that A c / <g; xq, so that 
n(x) « n(xo) throughout the duration of the interaction. 
Three phases of the interaction event are identified and 
illustrated in Fig. 2. 

1. Collision Phase (Fig. 2a). Both a forward shock 
(FS) and reverse shock (RS) are found in this 
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phase; the forward shock accelerates the cloud ma- 
terial, and the reverse shock decelerates the shell 
material. 

2. Penetration Phase (Fig. 2b). This phase bifur- 
cates into two cases, depending on whether the RS 
crosses the shell before the FS crosses the cloud, 
or whether the FS crosses the cloud before the RS 
crosses the shell. In the former case, the shocked 
fluid produced during the collision phase deceler- 
ates as more cloud material is swept up at the FS, 
during which time a new accelerated particle pop- 
ulation is introduced at the decelerating forward 
shock, so we denote this process the deceleration 
shock (DS). In the latter case, the shocked fluid 
produced during the collision phase is accelerated 
by the remaining shell material at the RS, and a 
new particle population is accelerated at the RS. 
We only treat the first case of this phase (which is 
the most important case for GRB studies) in this 
paper. 

3. Expansion Phase (Fig. 2c). After the FS crosses the 
remaining cloud material in the first case of Phase 
2, or the RS crosses the remaining shell material 
in the second case of Phase 2, the shocked fluid, 
being composed of highly relativistic particles and 
magnetic fields, expands and adiabatically cools. 

2.2.1. Collision Phase 



We follow the appro ach of ISari &: Piranl (|1995f) and 
iKobavashi et alj ([1999) to analyze this phase. Let T = 

1/yl — P 2 represent the relative Lorentz factor of un- 
shocked shell material in the rest frame of the shocked 
fluid. From Fig. 2b, 



f = rr (i - pp ) 



-(- 

2 v r 



in) 



where the last relation applies in the_ regime ro,r ^> 1 
considered here. Note that Y and Y remain constant 
during the duration of Phase 1, as a consequence of the 
assumption that A c ; <C x. 

The shocked fluid travels with Lorentz factor Y, and 
the shock itself moves with Lorentz factor Y s . When 
r>l, then r s = \/2Y, implying a compression ratio of 
4T (particle densities will_always be referred to the proper 
frame). Likewise, when r>l, the Lorentz factor of the 
RS in the comoving fluid frame is t r = y/2Y. When 
r - 1 « 1, then (3 r = y/\ - tr 2 = 4/3/3, implying a 
compression ratio = 4. Thus T r = \J~2T when f > \/2, 
and j3 r = 4,3/3 when f < \J~2~. The Lorentz factor_of the 
reverse shock in the stationary frame is r r = rr r (l — 

00 r ), and P r = y/l -Tr 2 . 

When r > 1, the fluid density of the FS material is 
nf = 4Tn c i . The density of the reverse-shocked fluid is 
n r = (4F + 3)n(x) for a relativistic reverse shock, and 
n r = (4 + hp 2 1 A)n(x) for a nonrelativistic reverse shock 
composed of cold shell material; an expression that joins 
the two regimes is n r rs 4Tn(x). The equality of kinetic- 
energy densities at the contact discontinuity implies that 




10" 10 10- 



10" 6 10' 4 



F/4r 



Fig. 3. — Lorentz factor F of the shocked fluid in the explosion 
frame, and Lorentz factor T of the unshocked shell material in 
the rest frame of the shocked fluid, as a function of F/4Tq when 
To = 300. Solid curves give numer ical resul ts, and dotted curves 
show approximations given in eqs. il 1 41) and l|15|l . 

The relativistic shock jump conditions fo r an isotropic 
explosion in a uniform CBM therefo re give (|Sari &: Piranl 
U995tlPanaitescu fc Meszaroslll999h 



i(xo) 



F 



En 



n cl ATTxlTln c im p c 2 A(x ) 

( 2T 2 /P 2 , Nonrelativistic RS (NRS) 
f 2 ^7 "* | t.2 /r2 



r 2 



T 2 /T 2 , Relativistic RS (RRS) 



Hence 



whenF> 4r^ (NRS), 



(13) 



(14) 



and 



v/2r - FV2 V 2 



when F < 4Fg (RRS) 



(15) 



Fig. 3 shows the shocked fluid Lorentz factor T and the 
relative Lorentz factor f as a function of F j 4Fq for Tq = 
300. The accurate numerical results are shown by the 
solid curves, obtained by solving F = P(T 2 — T)/P(T 2 — 
r). The dotted curves show the approximate analytic 
expressions, eqs. ([T4|) and (fl~5| . using the expression T = 



n/(r-l) =n r (r-l) S 4n c/ r 



4n(a;)(r 2 -r). (12) 



y / T Q /2F 1 / 4 in the latter case. The numerical solutions 
for r and T are used in subsequent calculations. 

The FS power is dE'/dt'\ FS = A cl n cl m p c 3 P(T 2 - 
T) = A c in c im p c 3 r 2 . The reverse shock power 
is dE'/dt'\ RS = A cl n(x a )m p c 3 P{t 2 - f). Thus 
(dE' /dt'\ R s)/(dE' /dt'\ FS ) = 1, so that equal power is 
dissipated as internal energy in the FS and the RS. The 
comoving duration of the FS before penetrating the cloud 
is Ai' FS = A c ; /Yc. The comoving duration of the RS be- 
fore penetrating the shell is At' RS w YqA/PYc. Thus 
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case (1), where the RS crosses the shell before the FS 
crosses the cloud, applies when At' RS /At' FS < 1, that 
is, when A/A c/ < (3T/T T < 1/T . Case (1) generally 
applies when A -C x/Tq (eq. [2]) and A c ; < x/Tq, which 
is the condition for STV. 

The range of x occupied by the FS fluid during phase 
1 (designated as FS1) at time t* = t z + xcosO/c is 

x + /3c(t* — t*o) < x < x + fi s c{t* — i*o) , (16) 

while the range occupied by the RS fluid during phase 1, 
or RSf , is 

x + p r c(t r — i* ) < x < x + [3c(t* — i* ) , (17) 



where £»o = xq/PqC and f3 s — vl— Ts • The comoving 
time t' — (t^o/To) + (i* — t*o)/T, so that the elapsed 
comoving time since the start of the interaction is At' = 
(i* — t*o)/T. In the absence of an interaction, the shell 
occupies the range poci* — A < x < poet*. Thus the RS 
leaves the shell when poci* — A = xq + pVc(i* — £»o), that 
is, when 

t» = <»A = t*0 + To oT~ ■ (18) 

(A) - Pr)C 

This represents the end of Stage 1 for the case where 
the RS crosses the shell before the FS crosses the cloud. 
Therefore Stage 1 lasts for comoving time 



A4s " Wo 



A 



(19) 



which may be compared with the approximation previ- 
ously derived. The time measured between the start of 
the GRB explosion and when the blast wave first encoun- 
ters the cloud is 



*o = (l + z)|°-x 



i-A) 



for 0; < 6 r 



1 - /3 cos(6» i - Od) for Oi > 6 cl . 

(20) 

Following the treatment of iSari et al.l (fl998l) (fo r the 
accuracy of this approach, see iDermer et alj 120001 ). we 
let n("f; x, t')d^ represent the differential number density 
of electrons with Lorentz factors between 7 and 7 + d'y in 
the shocked fluid. Assuming that all swept-up electrons 
are accelerated by the FS with spectral index pf, joint 
normalization of number and power gives a minimum 
electron Lorentz factor 



7mm,/ = e eJ ^-f(p f ) (r - 1) , where f(p) = (? 



p-V 
(21) 

and e e j represents the fraction of swept-up FS power 
carried by the injected nonthermal electrons (subscript 
"/" refers to the FS, and subscript V to the RS). For 
electrons accelerated by the RS with spectral index p r , 
similar considerations give a minimum Lorentz factor 



7n 



= £ e ,r — f(Pr) (r 



1 



(22) 



Electrons cool through synchrotron losses during co- 
moving time At' since the start of the collision to cooling 
Lorentz factor 



7c,/ 



67rm e c 
a T BjAt 



7 and 7 C 



67rm e c 
cr T B 2 r At' 



(23) 



for the FS and RS fluids, respectively, assuming domi- 
nant synchrotron losses. The magnetic field strength B 



in the FS and RS fluids is assigned according to the usual 
prescription, namely that the magnetic-field energy den- 
sity is proportional to the downstream energy density of 
the shocked fluid with proportionality constant eg < 1. 
Thus 

B 2 j = 327rn c im p c 2 eBjl3(r 2 — V) , and 

B 2 r = 327m(x)m p c 2 e Bir /3(f 2 - V) . (24) 

The maximum electron Lorentz factor jmax is deter- 
mined by equating the most rapid acceleration rate ex- 
pected in Fermi processes, 7 occ = tieBjm e c 1 where the 
rate factor £2^1, with the synchrotron loss rate |— 7 S j,n|, 
giving 



In 



,/ = ( 



67ree 2 ,/^i/ 2 ^ 1-2 x 10 8 e 2 ,/ 



(25) 



a T B f VBj{G) 

and similarly for jmax.r at the RS. 

We calculate the nonthermal synchrotron radiation in 
the collision phase spectrum by noting that the emissiv- 
ity in eq. ^ can be expressed through the relation 



e'j'(e',n';x,t') S -ca T U B j 3 s n(j s ;x,t') 



(26) 



where Ub = B 2 /8ir is the magnetic-field energy density, 

ls - VmW^) ' (27) 

and the critical magnetic field B cr = m 2 c 3 /eh = 4.414 x 
10 13 G. This follows from the formula for the electron 
energy-loss rate through synchrotron emission, and from 
the expression for the mean energy of synchrotron pho- 
tons radiated by an electron with Lorentz factor 7 S . Eq. 
(|26|) applies to the FS and RS fluids by using values of 
B and 7 S appropriate to each component. 

Assuming that the electrons are accelerated by the 
first-order Fermi shock process and injected into the 
shocked fluid in the form of a power law with index p, 
then the density distribution 71(7; x, t') of energized elec- 
trons can be approximated by the expression 

771(7; a;, i') = 
(7/7o) s_1 for7o<7<7i, 



n(x, t') 
s - 1 



(71/70) 1 s (7/7i) p , for 71 < 7 < 7„ 



(28) 

In the slow-cooling regime, 7 mi „ < j c , 70 = 7 mm , 71 = 
7c and s = p. In the fast-coolin g regime, 7 C < 7 mm , 
7o = 7c, 71 = 7mm, and s = 2 (jSari et al .l Il998h . At 
energies e < e = (B/B cr )jQdu/(^ + z), corresponding to 
7 < 7oj the synchrotron spectrum is approximated by the 
elementary synchrotron emissivity spectrum (oc e 1 / 3 in a 
vF v representation) . Thus we can include this branch of 
the synchrotron spectrum by writing 

Wt;^)^ (-) 1/3 for 7 < 70. 
s - 1 eo 

Synchrotron self- absorption is neglected here, which is 
valid for radiations detected at infrared frequencies and 
higher considered in this paper. 

The proper densities of electrons swept up by the FS 
and RS are 

f3Tn(x) 



n f (x,t') 



4rn c / and n r (x,t) = 



(P-Pr)r 2 ' 

(29) 
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respectively. In the treatment used here, the total elec- 
tron density is independent of location within the volume 
of the FS and RS fluids during Phase 1. The collision 
phase spectrum is calculated by substituting eq. (|26[) into 
cq. ([6]), making use of eq. ([27]) to relate 7 to the received 
photon energy. 

2.2.2. Penetration Phase 

Only the case where the RS passes through the shell 
before the FS passes through the cloud is analyzed here. 
During this phase, the shell is assumed to be hydrody- 
namically connected so that the entire shocked fluid de- 
celerates with the same Lorentz factor r 2 . As the FS 
passes through the remaining cloud material, it sweeps 
up and injects a new nonthermal particle population in 
addition to the nonthermal electrons left over from Phase 
1. We denote r = T\, the Lorentz factor of the shocked 
fluid during Phase 1. The equation for blast-wave de- 
celeration of a relativistic blast wave when it sweeps up 
material from the cloud is given by 

F2(X) = [1 + (2 - l p)T 2 1 c 2 m(x)/E ]y(^) (30) 

(jBottcher fc Dermedl2000t [D crmcr fc Humill200l . where 
ip is the fraction of swept-up energy that is radiated, and 
m(x) = Airrrip f x& dx' x' 2 n c i(x') = ^■nm p n c ix\(x — xa) 
is the swept-up mass. Here 

x A =x a + - — (31) 

PO — Pr 

is the distance of the FS from the explosion center at 
the end of Phase 1 (see eqs. [TB] and [TB]). Defin- 
ing the deceleration radius Xd through the expression 
Tim(xd) = Eq/Tic 2 gives the following expression for 
the deceleration radius for a blast-wave/cloud interac- 
tion: 

E 

Xd = xa + ^5 j- ■ (32) 

■iTrm p c 2 T(n c iX A 

Assuming adiabatic evolution (ip <C 1), eq. ([3H)l becomes 

r 2 = - Fl =- Fl (33) 

./l + 2(-2=2M y/l + 2p ls u-u 2 /2T 2 



where u = c(£» — t„ )/(xd — xa), and (3i s = y/l — 1/2T 2 . 
The expression on the right-hand-side of eq. ([33)) is ob- 
tained by integrating dx = /32s(x)dt*, where r 2;j (x) = 
V2T 2 (x). 

During Phase 2, the particles injected at the RS and FS 
during Phase 1 are no longer subject to further accelera- 
tion. We assume that these fluid volumes, denoted RS2 
and FS2, respectively, neither expand nor contract dur- 
ing Phase 2, and decelerate with Lorentz factor r 2 (i*). 
The range occupied by RS2 is therefore 

PrA , » PiA 

-+X 2 {t*-t*A) < X < X + - —+X2(t*-ti,A) 



x - 



00 ~Pr 



For FS2, the range is 
PiA 



(34) 



x - 



Po - P; 



+x 2 (i*-t*A) < x < XA+x 2 {t*-t*A) ■ (35) 



Here 



x 2 (t* 



(x d - x A ) 



,'■2 



r?-i 



-S ln| rYV^- 2 ^' + r? " 1 + S" 2/3lsl]| - 

(36) 

The elapsed comoving time since the start of Phase 2, 
obtained by integrating dt 1 = dt*/T2(t*), is given by 



XA r,U 

cTi U 2 



Pun 




PlsU 



2T 2 



p ls T 2 - v^rftarcsinGSi, - — j) - arcsin/3i s ]} . (37) 

The RS2 and FS2 emissivities are the same as the RSI 
and FS1 emissivities, except now At' = t' — t' A + (t»A ~ 
to)/T\ is used to evaluate the cooling Lorentz factor, 
eq. (|23|) . Moreover, because injection and subsequent 
particle acceleration has ceased in these two phases, the 
maximum electron Lorentz factor decays by synchrotron 
losses to a value of 



a T B 2 f 



(38) 



where jmaxj is given by eq. (|25p for FS2, with a related 
equation for RS2. Because the FS and RS fluid shells 
are assumed not to expand during Phase 2, the magnetic 
field remains the same as in Phase 1. 

The fluid containing the nonthermal electrons and pro- 
tons injected by the decelerating shock as it sweeps up 
cloud material in Phase 2, denoted DS2, occupies the 
range 



xa + x 2 (i* - **a) < x < xa + X2 S (t* 
where 



**a) , (39) 



V2Tl 



[ (|-A s r?) y - 2p ls u + 2T 2 - i+ / s la r? v /2rf - 1 ] 

(40) 

The DS2 emissivity at time t*, determined by the values 
of x, and observer time t z in the evaluation of eq. ([6]), is 
produced by nonthermal electrons injected at location xt 
and time t-n . The values of x, and t*j are determined by 
the intersection of x 2s = xa + X2 S (t*i — i*A) , represent- 
ing the leading edge of the deceleration shock (DS) that 
injects particles into the fluid, and X4 = x — x 2 {t* 
representing the worldline connecting the injection coor- 
dinates to the received values. The injection coordinates 
Xi and m = c(t*i — t*A)/{x d — xa) are therefore ob- 
tained by numerically solving x — xa — x 2 (i* — £*a) = 
X2s(t*i — <*a) — x 2 (i*i — £*a) using a staggered leapfrog 
routine. 

The received spectrum during Phase 2 is given by eq. 
©, with the x-integral limited to the ranges defined 
above. The emissivity for DS2 is given by eq. (|26p . with 
values of "f m in, 7c and "f m ax appropriate to the injection 
value of r(t*j). The elapsed comoving time since injec- 
tion used to calculate j c is obtained from eq. (j37|) . The 
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Doppler factor at time f* uses eq. (|33j) for T. Phase 2 
ends when x — X2, that is, at stationary time 



£*2 — i*A" 



\/2r x ( 



[ J2TI-1-J2TI - 1 - 2{- 2 



(41) 

2.2.3. Expansion Phase 

After the FS passes through the cloud, the shocked 
fluid has Lorentz factor 



r 3 = r 2 (t* 2 ) 



/l + 2(2a=2^) 



(42) 



Because there is no further pressure on the fluid, the 
shocked fluid is assumed to expand outward along the 
parallel direction with speed (3\\c and Lorentz factor 



1^ = 1/^1-/32. 

(Because the shocked fluid has a much smaller parallel 
than transverse extent, due to the compression of the 
quasi-spherical cloud material by a factor > 4r3 in the 
parallel direction, effects of transverse expansion can be 
neglected.) No further injection or acceleration takes 
place in Phase 3, and the relativistic electrons cool by 
synchrotron and expansion losses. The shocked fluid in 
Phase 3, denoted by RS3, FS3, and DS3, consists of fluid 
from the prior RS, FS, and DS phases. In order that 
each fluid phase is to have the same fractional volume 
expansion, the FS fluid is assumed to expand in both 
transverse directions in the comovin g frame with speed 
/3,| C /3, where ^13 = 1/^1-0311/3)*. 

From these assumptions, we define the four sta- 
tionary frame Lorentz factors Trs3 = || ( 1 — 

/3 3 /3||) , r FS3 - =r 3 r|| 3 (i-/?3/j||3) , r FS3+ = r 3 r p (i + 

/3 3 /3|| 3 ) and r_Ds 3 = r 3 F||(l + f3 3 f3\\) , along with their 
associated (3 factors (3rs3, Pf S3-, PFS3+ and PdS3, re- 
spectively, that characterize the motions of the radial 
boundaries of the different fluid layers. The range of 
RS3 at time t» is 



xo + 



pVA 

#>-A 

0i A 



+ :£2 (i*2 — t*A) + cf3RS3(t* — £*2) 



< x < x + +X 2 {t*2 - Ua) +C{3 FS 3~(U -U2) • 

PO - Pr 

(43) 



For FS3, the range is 
Pi A 



x 



+ X 2 (i*2 - i*A ) + c[3fS3- {U - t*2 ) 



<X<X A + X 2 (U 2 " **a) + c(3 F S3+(U - U 2 ) • (44) 
The range of DS3 is 

XA + £2(^*2 - Ua) + c(3fS3+(U - t*2) < X < 

xa + x 2s (U2 - Ua) + c/3ds3{U - U2) • (45) 
Let R'Jt') = R'P + 2/3\\ct' denote the comoving par- 
allel width of the entire shocked fluid, with additional 
subscripts i = "RS," "FS," and "DS" to refer to the 
widths R\Jt') = R',? t + §/?|| c(t' - t' 2 ) of the forward, re- 
verse, and deceleration shocked fluid layers, respectively. 



The superscript "0" refers to the radial width of the fluid 
layer at the end of Phase 2, and t' 2 is the comoving time 
at the end of Phase 2. From eqs. flS]) - (g5]), R'P RS = 

r 3 (p- x -pV)A/(/? -pV), R§p S = r 3 (j3 s -(3 r )A/(j3o-l3 r ), 
and Rf DS = T 



X2{t*2 — Ua)], with RiP 



3[X2s(t*2 - Ua) 
equal to the sum of these three widths. 

Conservation of magnetic flux for the transverse 
magnetic-field component B^ implies B±R' 2 = const — 

B Q ± Rf 2 , so that B x (t') = B Q ± Rf 2 /[R\P + |/3||c(t' -t' 2 )} 2 . 

If the electrons are isotropized on time scales short com- 
pared to the cooling timescale, then the electron energy 
loss-rate due to adiabatic and synchrotron losses is given 

by 



G?7 



1 dR' H cT T [Bl(t') 



7 + 



Bl{t')] 



6irm e c 



(46) 



where B\\ is the parallel magnetic-field component. As- 
suming that B±(t') ^> Bu(t'), which would occur if the 
magnetic field is formed by sweeping up and compressing 
an external field, then eq. (|46|) can be written in the form 



(47) 



dl 7 . ,7 2 
- = - + b~, 



dr t t" 1 

where r = 1 + |/3||c(t' - t 2 )/R\? and b = 

R'°o-TB° ± 2 /(4 Tr(3i\m e c 2 ). This can be solved 

(|Gupta et al.1 [2006) using the substitution y = "fT~ n , 
with n = 3, giving the result 



7(r) 



4r 3 



6(r 4 - 1) + 4r 4 / 7l 



(48) 



where ji is the initial electron Lorentz factor at t' = t' 2 
or t = 1. 

The emissivity for RS3 and FS3 is given by eq. (|26p 
with 



n(r,t') = n(r, t) = 



^11(0, -1^(7*5 1K 2 _ «(7,;i)7 2 



R 



1 



r 2 7 2 



(49) 

B —>■ B±(t') = B\lr 2 '', and n( 7 i; 1) is the electron 7 dis- 
tribution at the end of Phase 2. This expression also 
applies to the DS3 emissivity, with the emission coor- 
dinates x, i* mapped back to the injection coordinates 
Xi,t*i as was done for DS2. Note the minor inconsis- 
tency for Phase 3 between the original assumption that 
the magnetic field is randomly ordered, compared to the 
dominant perpendicular magnetic field component im- 
plied by field compression in this phase. This inconsis- 
tency should be relaxed with the use of directional syn- 
chrotron emissivity spectra in more detailed treatments 
but, in either case, the contribution from the expansion 
phase is usually much smaller than that from the pen- 
etration phase. The Doppler factor during Phase 3 is 

Su = [r 3 (i - /M]" 1 - 

3. RESULTS 

Nonthermal synchrotron emission spectra from blast- 
wave/cloud interactions were computed using the formu- 
lae given in the previous section. Standard parameters 
are given in Table 1. For the model GRB pulse, we 
consider a source at z = 1, corresponding to the mean 
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plied by the idealized cloud location and geometry. From 
eq. (|20|) for on-axis clouds, 



Fig. 4. — Lightcurves calculated at 511 keV for a standard model 
GRB pulse with eg = 10 — 4 and r] = l/To- The cloud geometry is 
fixed with 9 c i = 3/To, but its density is varied as labeled by the 
curves. Top figure: linear scale. Bottom figure: logarithmic scale. 



rcdshift of the pre-Swift BATSE/Beppo-SAX/HETE- 
2/INTEGRAL sample. The coasting Lorentz factor and 
apparent isotropic energy release of the blast wave are 
chosen to be T = 300 and E — 10 53 ergs, respectively. 
The initial width of the shell is Ao = 10 7 cm. The blast 
wave is assumed to encounter a small cloud with nominal 
density n c i — 10 3 cm -3 located at a distance xo = 10 16 
cm from the explosion center 3 centered along the line- 
of-sight to the observer. To simulate a quasi-spherical 
cloud, the cloud width A c / = 2 x 10 14 cm and 6 c i = 0.01. 
Other standard parameters are injection electron index 
p = 2.5, e e = 0.1, and ti = 0.1. The parameters e e ,e_B, 
£2, and p of the FS and RS are assumed to be the same 
throughout all phases of the interaction in the calcula- 
tions presented here. 

Fig. 4 shows lightcurves computed at 511 keV for the 
standard model GRB pulse with eb = 10~ 4 . We prove 
in Section 4.2 that a rapidly variable bright pulse cannot 
be made when r\ = 1, and numerical simulations with 
r] = 1 confirm this result. Thus we adopt the thin-shell 
assumption (r] = 1/ro) in Fig. 4, and assume that the 
cloud is on the axis to the line-of-sight to the observer. 
This calculation shows the dependence of the behavior 
of the light curve on cloud density n c i. 

The kinematic minimum and maximum times are im- 

3 The total mass of such clouds distributed in all directions 
around the GRB source is < 10 — Afm. 



tr. 



_ (1 + z)x 



2Tlc 



= 3.70 



1 + z 



r 2 

1 300 



s and 



£mm(l + Tq^;) 



(50) 

here T3oo = Fo/300. Depending on observing energy, 
the maximum emission time will also be extended by the 
timescale for electron cooling. As can be seen from Fig. 
4, at low densities the flux is weak and the blast wave 
hardly decelerates. For most purposes, such interactions 
can be neglected as they extract only a small amount of 
energy from that solid angle element of the blast wave. 
With increasing n c i , the vF v flux increases until the cloud 
density becomes so thick that the RS becomes relativis- 
tic, r <§C To, and the emission is received at lower photon 
energies. At such high densities, the flux received by a 
distant observer is weak but may last for a long time, 
especially at small photon energies. 

The model light curves show a sharp rise followed by 
a more leisurely power-law temporal decay. Some of 
this rapid rise has to do with the idealized cloud ge- 
ometry: the use of an actual spherical cloud geometry 
(which could be mimicked by superposition of interact- 
ing annuli) would tend to reduce the sharpness in the 
rising phase of the GRB light curve. As can be seen 
from the bottom figure of Fig. 4 in a logarithmic rep- 
resentation, the power-law decay of the flux displays a 
temporal curvature or softening at times well past the 
peak from cooling breaks in the synchrotron spectrum. 
This is related to the detection of higher energy pho- 
tons that were emitted from a higher-energy and usu- 
ally s ofter part of the sp e ctrum by off-axis emittin; 
sites (Fcnimorc ct al. 1996; Kumar & Panaitescu 200 
iDerme r 2004; Liang et a l. 2006; Zha ng et al.ll2007f ). The 



pulse shapes with n c i <C 10 cm cor respond to kine- 
matic "curvature" pulses |Dcrmer 2004) , where the radi- 
ating shell is thin and the emission is radiated promptly 
on the crossing time defined by the comoving shell thick- 
ness, and the spectrum is well described by a single power 
law in the weak cooling regime. The GRB light curve 
produced by interactions with such clouds with optimal 
density n c i = 10 3 - 10 5 cm -3 would represent a generic 
long-duration GRB except that it is rather weak, reach- 



2 x 10 8 ergs cm 2 s 1 



at 



ing peak vF v fluxes of ~ 1 
~ 500 keV. 4 

The weakness of the flux is largely due to the small 
value used for the es parameter. Fig. 5 shows the vF v 
SEDs calculated near the peak of the e = 1 light curve 
at observer time t = 4.46 s, for three different cases of 
€b and with n c \ = 10 3 cm~ 3 . Other parameters of the 
event are the same as in Fig. 4. The various components 
making up the spectrum are shown. The DS compo- 
nent during Phase 2 makes the dominant contribution 
to the vF v spectra near peak energy release. For the 
higher magnetic field case es — 0.01, and even more so 
for €b — 0.1, the stronger cooling causes the Collision 
Phase components to decay away more rapidly. 

For the case es = 10~ 4 on the left, one sees that the 
peak vF v flux, which is important for triggering, is at the 

4 The uF v peak flux measured by a GRB telescope integrating 
over a finite bandwidth would be a factor of a few larger due to a 
bolometric spectral correction. 
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Fig. 5. — Separate spectral components and the total vF v flux 
observed at 4.46 s after the start of the GRB, near the peak flux 
of its lightcurve. This GRB pulse model uses parameters given 
in Table 1, with ra c ; = 10 3 cm -3 and r) = l/To- The differences 
between the system are that, from left to right, eg = 10 -4 , eg = 
10 — 2 , and eg = 0.1. The total spectrum is given by the heavy 
solid curve with data points, the FS components are given by the 
solid curves, the RS components by the dashed curves, and the DS 
components by the dotted curves. The heavy blue, medium red, 
and light green curves correspond to Phases 1, 2 and 3, respectively. 

level of a few x 10 -8 ergs cm -2 s _1 . BATSE triggered 
between e = 0.1 - 0.6 at a typical vF v flux level of w 
5 x 10 -8 - 10 -7 ergs cm -2 s" 1 , so such a GRB would 
be hard to detect unless it occurred nearby (z < 1) and 
happenened to be pointed at us. This cannot be excluded 
for GRBs without rcdshift measurements, so our interest 
for spectral modeling is principally on GRBs that have 
rcdshift information. Stronger magnetic fields produce 
brighter pulses, and when > 0.01, the model pulse 
flux is at the level > 3 x 10 -7 ergs cm' 2 s _1 , making 
it an easily detectable GRB with BATSE at z = 1, and 
even more so when es ^ 0.1 and e e > 0.1, although 
the adiabatic assumption for the deceleration shock may 
then become invalid, a situation that can also arise even 
when t e . y 0-1 if photohadronic energy losses are strong 
(|Dermerll2007t ). 

An interesting result of these calculations is that the 
spectral index below the vF v peak is very hard, in num- 
ber index harder than —0.7, — 1, and —4/3 for es — 
10~ 4 ,10 -2 , and 10 _1 , respectively. In no case does a 
well-defined cooling spectrum with number index —3/2 
or softer form. This suggests ano t her a rgument against 
the concerns of Ghis ellini et alj (|2000f ) that radiative 
cooling spectra should be observed in GRB spectra. Here 
there is not sufficient time until the interaction is over for 
the electrons to cool, so the spectra remain much harder 
than a cooling synchrotron spectrum below E p k- What 
stops the interaction and cooling is simply that, after 
both shocks have passed through the shell and cloud, 
the disturbed material rapidly expands, shutting off the 
subsequent synchrotron emissions as nonradiative adia- 
batic losses increase and the magnetic field intensity, and 
therefore the synchrotron losses, decrease. 

The spectra formed by the interaction of a relativis- 
tic blast wave with a small cloud consisting of both FS 
and RS emission have the characteristic "Band" -shape 
(jBand et alj I1993D . except for a lower energy emission 
component from the RSs. The relative importance of 
the RS and FS components increases with increasing eg 
for systems with nonrelativistic or weakly relativistic RS. 
When eg S> 10 -4 , the RS produces a strong optical flux. 
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Fig. 6. — Light curves at e = 1 (hu = 511 keV) from the 
interaction of a GRB blast wave with a cloud for standard GRB 
pulse model with eg = 10~ 2 , r\ = 1/To, and other parameters 
given in Table 1, for different values of n c i. Top figure: linear 
scale. Bottom figure: logarithmic scale. 



Indeed, this is the external-shock model explanation for 
prompt optical emiss i on in GRB emission, first discov- 
ered bv lAkerlof et all (| 1999( 1. 

Even when cb — 10~ 4 , the RS emission produces con- 
cavity in the X-ray spectrum that could be detectable in 
joint spectral fitting of the Swift BAT and XRT data with 
optical/UV measurements with the UVOT and ground- 
based robotic optical telescopes. For GRBs as modeled 
by Fig. 5, optical flashes as bright as my w 10 - 17 would 
be coincident with GRB pulses. The relative amplitude 
of the reverse shock component is also, however, very 
sensitive to the value of the shock thickness parameter 77, 
and various underlying assumptions for the constancy of 
the e e and €b parameters with time and equality in the 
reverse and forward shocks. The rich parameter space 
could allow detailed modeling of prompt and early after- 
glow optical light curves from, e.g., RAPTOR, ROTSE, 
and super-LOTIS. 

Fig. 6 shows 511 keV light curves calculated for the 
standard GRB pulse parameter set with €b = 10~ 2 , 
while varying the cloud density n c i. The same basic be- 
havior as was found in Fig. 4 is apparent, except that the 
pulse is much brighter and decays more rapidly. The op- 
timal cloud density to make the brightest possible pulse 
is again near n c i = 10 3 cm -3 for €b = 10~ 2 (see Section 
4.2). The generic pulse profile found here is typical of 
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Fig. 7. — Calculated light curves at e = 1 (hv = 511 keV) 
from a GRB blast wave impacting a cloud with opening half-angle 
9 c i = 0.01 located on the symmetry axis of the jet and the line- 
of-sight to the observer. Parameters for the standard GRB pulse 
model given in Table 1 are used, with eg = 0.1 and r) = 1/To. The 
extent of the jet opening angle is 8j = 0.01,0.003, and 0.001, as 
labeled. Top figure: linear scale. Bottom figure: logarithmic scale. 
The straight line is oc t~ 1 . 

a "FRED" -type GRB. 5 Such a profile could just as eas- 
ily be produced by a GRB bl ast wave deceleratin g in a 
uniform surrounding medium (|Dermer et al. 1999). The 
external shock model can explain relatively smooth light 
curves that moreover show various widths and asymme- 
tries. The central question for GRB pulse modeling is, 
however, whether the rapid variability found in some 
GRB light curves can be reproduced in an external shock 
model. 

Fig. 7 displays the light curve formed for the standard 
parameter set with eb = 10 _1 and r\ = 1/ro. Here 
the cloud has extent c i = 0.01, and a jet with opening 
half-angle 9j collides with the cloud. This plot indicates 
the latitudes from which the pulse flux originates. The 
high-latitude emission arrives at late times, and most of 
the fluence is formed by interactions for jet interactions 
within the Doppler opening angle 1/To. Although this 
calculation illustrates the part of the blast- wave/cloud 
interaction from where the flux originates, the opposite 
limit, where 0j 3> 9 c i, is of more interest for producing 
STV in GRB light curves. 

Fig. 8 shows the emission profiles produced by clouds 
with density n c i = 10 4 cm -3 and typical sizes A c / = 

5 FRED stands for Fast Rise, Exponential Decay, though the 
decay phase is often better fit by a power law. 
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Fig. 8. — Pulses formed by GRB blast wave shells with r\ = 
1/ro colliding with clouds at different inclinations from the line of 
sight. The clouds have half-opening angle 8 c i = 0.001 which, for 
To = 300, represents an angular extent equal to 1/3 of the Doppler 
cone. Top figure: linear scale. Inset shows detail of brightest pulse. 
Bottom figure: logarithmic scale. 

0.001 xio = 0.3xo /Tq = 10 13 cm, so that their angular 
extent is smaller than the Doppler angle. The blast wave 
is assumed to form a thin shell with r\ = 1/ro. I n an ide- 
alized circumburst medium where cloud properties were 
roughly the same throughout this region, then for every 
bright pulse or two, there would be another ~ 7 dimmer 
by a factor of 10 - 20, and another 20 pulses dimmer than 
the brightest pulse by a factor ss 100. Pulse widths of a 
second or so are formed by the brightest p ulses for such 
cloud s, in keeping with the pulse paradigm (|Norris et al.l 
1996). At late times and from larger angles a lower flu- 
ence emission plateau would be formed if the jet and 
clouds are assumed to maintain roughly uniform proper- 
ties at these angles. 

Because of the unknown distribution of cloud densi- 
ties, sizes and locations, as well as effects of GRB jet 
structure, simulating GRB light curves within an exter- 
nal shock scenario is very model-dependent. But simu- 
lations of GRB light curves with small clouds randomly 
distributed in a she ll show many of the featu r es of spiky 
GRB light curves (|Dermer fc Mitmanl 119991 12004 and 
below, Section 4.5). 

The calculations shown here mean that given the 
thin-shell approximation i] = 1/ro, essentially all 
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Fig. 9. — Light curves calculated at e = 0.01 [hv as 5 keV) 
from a blast- wave shell in the thin-shell approximation r] = 1/To 
for a blast-wave Lorentz factor To = 100 impacting an on-axis 
stationary cloud with opening angle 6 c i = 0.01. Top figure: linear 
flare profiles as a function of n c ;. Inset: SED measured at observer 
time of 512 s with n c i = 10 3 cm -3 . Bottom figure: logarithmic 
flare profiles as a function of n c i . 

BATSE GRBs without redshift information, and many 
BATSE/Beppo-SAX/HETE-2 GRBs with redshift infor- 
mation can be modeled within an external shock sce- 
nario. The standard jet properties are Eq = 10 53 
ergs and z = 1, in accord with pr e-Swift observations 
of GR Bs with redshift information ([Friedman fc Blooml 
|2005| ). These pulses h ave SEDs in agre ement with 
BATSE observations (e.g. lPreece et alj |2000). bearing m 
mind different ways to r emedy the "line-of-death" prob- 
lem (jPreece et all 11998( 1 in a nonthermal sy nchrotron- 
shock model through the jit ter mechanism (iMedvedevI 
2000), ra diation reprocessing [iDermer fe Bottcherl l2000) , 
etc. (see IZhang fc Meszaroa 120021 . for a discussion of 
GRB spectral break models). 

Many X-ray flares observed with Swift can also be 
modeled within the thin-shell approximation, as illus- 
trated by the calculation shown in Fig. 9. Here the stan- 
dard X-ray flare parameters from Table 1 are used. The 
principal difference between these parameters, which ap- 
ply to Swift observations, and those for the BATSE GRB 
pulse observations is that the source redshift is z = 2 
and the apparent isotropic energy release is 10 ergs, 
as implied by o bservations of GRBs at such redshifts 
(|Ghirlandall200l . The flares were detected by the XRT 
on Swift, so that the mean telescope photon energy is 
e < 0.01. As can be seen, flares at flux levels reach- 
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Fig. 10. — Light curves calculated at e = 0.01 (hv i=s 5 keV) from 
a blast- wave shell/stationary cloud interaction in the frozen-pulse 
approximation r\ — > for a blast-wave Lorentz factor To = 300 
impacting an on-axis cloud with opening angle 9 c t = 1/300 and 
Aq = 10 9 cm. Top figure: linear flare profiles as a function of 



incidence angle 0i for n c ; 
flare profiles. 
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Bottom figure: logarithmic 



ing 10 -9 ergs cm~ 2 s _1 can be made by external shocks. 
Comparison with observatio ns of GRBs with me asured 
redshifts in the Swift catalog (jO'Brien et al. 2006) shows 
that most flares are well below this level, except for a flare 
from GRB 050820A at z = 2.612 that reached a vF v flux 
between 0.3 and 10 keV of w 10~ 8 ergs cm -2 between 
200 and 300 s after the GRB. 

The SED for the model GRB X-ray flare peaking hun- 
dreds of seconds after the start of the GRB is shown by 
the inset in Fig. 9. In these flares the reverse-shock syn- 
chrotron component can make an enhanced optical/IR 
emission component due to the concavity produced by 
the various RS components. 

GRBs which lack redshift information, such as GRB 
050502B with its enormous X- ray flare exceeding 10~ 9 
ergs cm -2 ( Falcone et al . 2006a), pose no difficulty to the 
external shock model. 6 If the frozen pulse assumption 
is allowed, then there is no difficulty in explaining 7- 
ray pulses or X-ray flares in GRBs, provided that the 
surrounding medium of a GRB is clumpy on well-defined 
density and size scales. 

A straightforward prediction of the external shock 

6 Because the vF v peak flux was above the Swift BAT energy 
window, the fraction of total GRB fluence in the X-ray flare is 
unknown. 
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model as presented here is that the spectra of X-ray flares 
tend to become softer with time, and will be systemat- 
ically softer than the pulses in the prompt phase. This 
is because the single expanding blast wave becomes less 
dense with radius, so that interactions with surrounding 
cloud material tend to have weaker forward shocks. Thus 
there will be a systematic reduction of the vF v peak ener- 
gies and a softening of the X-ray flare spectra. Softening 
of spectra are obse rved in flares from i ndividual GRBs 
like GRB 050502B (|^Icone_glaLj[2006i) or GRB 050822 
(jGodet et al.ll2007t ). This cannot be considered a defini- 
tive test of the external shock model, however, because 
such behavior can also be accounted for by a colliding 
shell/internal shock model if the relativistic ejecta tend 
to have lower T factors or less energy with time. It would 
be surprising and contrary to expectations of this model 
if X-ray flares tend to be harder than the prompt emis- 



4. PHYSICS OF PULSE AND FLARE FORMATION IN THE 
EXTERNAL SHOCK MODEL 

We find from this analysis that short pulses or flares 
can be made if the GRB blast wave undergoes no trans- 
verse spreading and interacts with small dense (~ 10 3 
cm~ 3 ) clouds within w 10 16 - 10 17 cm from the GRB ex- 
plosion. The formation of strong forward shocks with 
r = r is necessary to make bright pulses, and was 
the underl ying assumption of t he GRB lig ht curve sim- 
ulations in lDermer & Mitmanl (|1999l . l2004h . where GRB 
light curves similar to actual GRB light curves were sim- 
ulated with ~ 10% efficiency The implications and plau- 
sibility of the assumption of a thin blast wave shell are 
now considered. 

4.1. Temporal Behavior 

The characteristic timescales shown in the pre- 
ceding figur es can be understoo d from simple con- 
siderations (iFenimore et al.1 119961: ISari k, Piranl 119971 : 
iDermer fc MitmaiJll999t lloka et al.ll2005f) . Eg. (|9| gives 
the time 

[ ^ r ' )X {l-M, (51) 



t 



foe 



since the start of the GRB explosion when emission is 
first received from a cloud at location x and inclination 
angle Qi = arccos that is impacted by a blast wave 
moving with Lorentz factor r . The angular timescale 
giving the temporal duration over which emission, emit- 
ted at radius x in the angular range 02 = 8% + 6 c l > $i — 
max(0, 9i — 9 c i), is received is simply 



when 



At ang x 
1 + z foe 

: > 6 c i, and 

At ang 



(Mi - M2) 

\-e\) ~ 



2<i, r >i 



1 + z 2c 



(52) 



(53) 



when 9i < 6 cl . 

The radial timescale for emission from two points at 
different distances x\ and 22 but at the same inclination 
angle 8i is given from eq. (j20f by 

A cl 



where t m i n is defined by eq. (|50p , and the inclination and 
cloud angles are in units of the Doppler angle 6>o = 1/ro 
by the relations 

Vj = T 9i and v c i = 9 c iT , 

respectively. The characteristic duration of a shell/cloud 
collision is therefore 

2v c / . n\ 

At = tmin max[-j^(l + v i ), 

4viV c /i7(v 4 - v c i) + fa + v cl ) 2 H(v c i ~ Vi)] , (55) 

where H{u) = 1 if u > 1 and H (u) = otherwise. Note 
that the angular timescale is independent of cloud loca- 
tion x, using v c ; = TqOci = (A c ;/x)ro in eq. (|55|) . but 
depends on cloud size A c / and inclination angle The 
radial timescale is also independent of x when written in 
terms of cloud size. 

The STV condition for the external shock model is 
v c / <^ 1. Thus we see that the kinematic duration 
is determined by the angular timescale everywhere ex- 
cept when Vi < l/2Fo, that is, for nearly on-axis clouds 
with 9i < 1/2Fq. The kinematic duration of the pulse 
At oc 1/0; for 6i > 1/2TI. The pulse shortening for 
nearly on-axis events, combined with beaming factors 
and spectral effects, is shown from analyses and Monte 
Carlo simulations capable of expla ining the range of 
observed G RB BATSE light curves (jDermer fc Mitmanl 
fl999l[200l . 

4.2. Optimal Cloud Density for Bright Pulses 

The optimal density giving the brightest measured flux 
can be derived from three requirements: 



1. A strong FS which, from eq. (|13| . implies a maxi- 
mum cloud density given by 

Eq 



n c i 



< 



16tt x 2 TQm p c 2 A(x) 



(56) 



2. Significant blast- wave deceleration to provide effi- 
cient energy ext raction, which requires cl ouds with 
thick columns (|Dermer fc Mitmanl fl999). that is, 
with densities 



n cl 



> 



En 



4irx 2 r^m p c 2 A c i 



(57) 



3. The requirement of a strong FS and a thick column 
therefore translates into the requirement that 

A cl > AT 2 A(x) (58) 

in order to produce STV. The STV condition is 

(59) 



_ Acl_ _ Vci 1 

x r r 



(54) 



Using eq. © for the shell width, eq. ([55)) becomes 

4^0 < + 4 V T « 1 , 

x 

so that the requirement on the width of the radiating 
shell material in the external shock model for rapid vari- 
ability is simply that 

V Z 4r when x > 4r 3 A S 10 15 r 3 00 A 7 cm . (60) 

41 o 
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It is therefore not possible to have STV if 77 = 1, be- 
cause then the reverse shock is relativistic, the forward 
shock Lorentz factor is low and Doppler boosting is lost. 
But it becomes possible in the thin-shell approximation 
77 ~ l/IV As shown by the light curves and SED calcu- 
lations, the phenomenology of GRBs can be explained in 
the frozen pulse (77 = 0) and also, possibly, the thin-shell 
approximation. 

The optimal cloud density for bright variable pulses 
in the external shock model, from eq. (|56p for a GRB 
blast wave with apparent isotropic energy release Eq = 
10 54 E 54 ergs and shell width A{x) = 10 7 A 7 cm that 
encounters a cloud at 10 16 xi6 cm from the center of the 
GRB explosion, is 



n c i(cm~ 6 ) « 1.6 x 10 b 



E< 



(61) 



2 p2 A 

For instance, the calculations shown in Fig. 5 with r/ = 
1/To have A 7 (a; = 10 16 cm) = 38. 

4.3. Width of Radiating Shell 

Besides the requirement for STV that the unshocked 
fluid shell be thin or frozen, the width of the shocked fluid 
shell cannot be too great, for otherwise rapid variabil- 
ity would be washed-out by contributions from different 
parts of the shell. 

First consider a blast wave sweeping up matter from a 
uniform surrounding circumburst medium with density 
no(cm~ 3 ). The swept-up material, compressed by the 
factor 4r at the external shock, has thickness 

Khto = yff and A sh (x) = ^ (62) 

in the proper and explosion frames, respectively. Because 
the angular extent of the shocked fluid layer in the blast 
wave covers the full Doppler cone, the variability is al- 
ready limited to a timescale > 0.25i m i„ by light travel 
time limitations. Under these circumstances, rapid vari- 
ability is not expected, even when considering density 
jumps (TNakar k Granotll2"00l . 

Now consider the transverse width of the shocked 
fluid shell colliding with a cloud with Dopper size 
v c i <C 1 that subtends the solid angle Af2 c ; as 
seen from the center of the explosion. By equating 
the mass m p AQ c ix 2 n(x)A' sh (x) with the cloud mass 
m p AQ c ix 2 n c iA c i we have, using the relativistic density 
jump condition n(x) = 4rno, 

A, 



= = ^ and A s ^) = ^ , (63) 

i.e., a shocked-fluid layer thinner than required to pro- 
duce STV. 

4.4. Thin Shell Assumption 

From the preceding discussion, one can see that the 
viability of an external shock model for the 7-ray pulses 
and X-ray flares depends on whether the GRB blast- 
wave width spreads in the coasting phase according to 
cq. ( |2l), with 77 < 1/Tq- In the gas-dynamical study 
of Mcszar os et al.l (|1993f ) . inhomogeneities in the GRB 
fireball produce a spread in particle velocities of order 
\v - c\/c ~ T~ 2 , so that A(x) ~ x/T^ when x > T^Aq 
and r\ ~ 1. This expressi on was also obtain ed in the 
hydrodynamical analysis by IPiran" et~ai1 (fl99l . 



Several points can be made about these results. First, 
gas-dynamical or hydrodynamic analyses omit the fluid 
nature of the explosion plasma whose behavior depends 
crucially on the magnetic field and MHD properties of 
the fluid. Second, the spread in A considered for a spher- 
ical fireball is averaged over all directions. As the fireball 
expands and becomes transparent, the variation in fluid 
motions or gas particle directions over a small solid angle 
~ 1/Tq of the full sky becomes substantially less. Small 
values of r\ could be realized if there is little dispersion 
in the baryon-loading of the fluid shell within the small 
range of solid angles subtended by the jet from which 
the observed emission is radiated. Third, the particles 
within the cold unshocked blast-wave shell will expand 
and adiabatically cool so that the fluid will spread with 
thermal speed Vth — PthC- The comoving width of the 
blast wave is r A + f3 t hcAt' m T A + p t hx/T , so that 
the spreading radius x spr = TqAq/ f3th- Adiabatic expan- 
sion of nonrelativistic particles can produce a very cold 
shell with (3q < 10 -3 , leading to very small shell widths. 

The requirement on the thinness of A(x) does not ap- 
ply to the adiabatic self-similar phase, where the width of 
the shocked fluid shell is A s h ~ x/12 Fn, as implied by the 
relati vistic shock jump equations (|Blandford & McKed 
[197(1 : see eq. (|62p . Consequently, no extreme variability 
is expected when a blast wave, already in the adiabatic 
self-similar regime, encounter s a den sity jump, as demon- 
strated by iNakar fc Granot] (|2007l ) ( in agreement with 
the results of unpublis hed numerical simulation s using 
the code described bv IChiang &; Dermerl (|1999l )). But 
even in this case, however, A <C x/T 2 if the blast wave is 
highly radiative (|Cohen et al.lll99a ). which could result 
from leptonic processes when e e ,es > 0.1, or hadronic 
processes when e e ~ 0.1 and cb > 0.1 and the surround- 
ing circumbur st medium is dense, with no > 10 3 cm -3 
(|Dermerll20d7r i. 

This question deserves considerably more study, and 
will require at least an MHD model employing self- 
consistent particle distributions and magnetic field ge- 
ometries in the erupting plasma, generation and amplifi- 
cation of the magnetic field in the shocked fluid layer (in- 
cluding particle acceleration and diffusion into the cold 
fluid shell), and magnetic field and particle energy evo- 
lution from shell expansion. 

4.5. Density Contrast 

A further requirement on this model is that the dense 
clouds are found in a medium sufficiently tenuous that 
the blast wave has not undergone significant deceler- 
ation by sweeping up this material. This constraint 
can be expressed by the condition that the deceleration 
radius Xd = (iE^/Anrn. D c 2 n ^T 2 ) 1 ^ ([Rees fc Meszarosl 
119921 iMeszaros fc Reeslll993| ) be much greater than the 
radius defining the condition for a thick column, eq. (|57|) . 
This translates to the requirement that the density con- 
trast between the surrounding medium density no and 
the cloud density n c ; is given by 



no 



n c i V E /mpC 2 



i(r 10 r 300 



> c; /10 3 cm 3 )(A c/ /10 9 cm) 3 



E 



54 



(64) 
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Fig. 11. — Simulated GRB lightcurves from the interaction of a 
thin blast wave with clouds with thick columns, using parameters 
from Table 2. The top four figures show ^ 100 keV light curves. 
The lower curve is the GRB emission and the upper curve is this 
emission added to Gaussian noise typical of the BATSE detectors. 
The bottom two figures show light curves at 5 keV characteristic 
of X-ray emission detected with the XRT on Swift. 

This is a severe constraint that can only be realized in an 
extreme environment such as the one described in more 
detail below. 

4.6. Simulations of GRB Light Curves 

Fig. 11 sho ws Monte Carlo simulation s, described in 
more detail in iDermer fe Mitmanl ((2004), of GRB light 
curves under the assumption that the cold blast wave 
shell remains thin and the interaction forms a strong for- 
ward shock. The parameters of the simulation are shown 
in Table 2, and here we use fixed cloud sizes with radius 
r\ that are "uniformly randomly" distributed in the vol- 
ume between R\ and R^. In all cases, the surface filling 
factor £ = 10%, and a fraction ( = 10% of the blast 
wave energy intercepted by the cloud is transformed into 
radiation received by the detector at energy m e c 2 e. The 
angular extent of the jet is kg/To, and an occulting fac- 
tor is included so that the portion of the blast wave that 
intercepts a cloud at a smaller radius no longer radiates 
when interacting with clouds at larger radii. 

The top four figures represent the GRB prompt phase 
detected near 100 keV using parameters appropriate to 
BATSE, including noise at a level typical of the BATSE 
detectors. The bottom two figures make use of parame- 
ters that would apply to Swift. The variety of simulated 
light curves is endless, depending on the seed used in the 
random number generator and the choices of cloud and 
blast wave parameters. Precursor events can be made if 
some material is found very close to the GRB (|Piro et al.l 
2005, compare figure E). 



4.7. Criticisms of the External Shock Model 

We describe the general criticisms that the external 
shock model for prom pt GRB emissions has endured. 
iFenimore et al.1 (|1996t ) demonstrated that an expanding 
blast wave would, if briefly illuminated over a large por- 
tion of its surface covering the Doppler cone, make pro- 
gressively longer pulses in accordance with the relation 
t = (l + z)x/T 2 c. The tendency for pulse lengthening can 
be avoided if a small fraction of the blast wave surface 
area within the Doppler cone were illuminated at any one 
times, but this was thought to require an unreasonable 
number of illuminating regions ( "clouds" ) . 

One approach is to estimate a maximum "surface fill- 
ing factor" of radiating sites consistent with measured 
varia bi lity, using kinematic arguments (jFenimore et al.l 
11999ft - ISari fc Piranl (|1997l ) argued that highly variable 
light curves cannot be constructed from blast-wave in- 
teractions with surrounding density inhomogeneities by 
superposing pulses made near the outer region of the 
Doppler cone to show how they overlap due to the large 
angular timescales; hence the sources of GRBs must re- 
quire an active cen t ral en gine to make discrete pulses. 
IDermer fc Mitmanl (|2004h showed that this argument, 
made also more recently by iPiranl (|2005h and lloka et al.1 
(2005), mistakenly uses the angular time scale near the 
outer edge of the Doppler cone, and misses a number 
of important points: One is that the beaming factors, 
not considered in these papers, heavily weight the total 
fluence for near on- axis clouds, and even more so the 
pulse peak flux because of the smaller kinematic vari- 
ability timescale for nearly on-axis clouds; another is the 
different (observer) times for the clouds nearly on-axis 
as compared with those off-axis, so that flux ratios from 
different parts of the shell are time-dependent; a third 
is the assumption that clouds are distributed with az- 
imuthal symmetry within the Doppler cone; a fourth is 
that the beaming factor of the relativistic jet that forms 
7-ray pulses in the prompt phase may be smaller than 
~ r^ 1 and the typical angles inferred from optical af- 
terglow breaks, which will reduce the importance of off- 
axis events. The crucial role of the narrow cold blast- 
wa ye fluid shell (no t the s hocked fluid shell considered 
bv iNakar fc Granotl |2003)) to ensure a strong forward 
s hock has also been overlook e d in p ast studies. 

iRamirez-Ruiz fc Fenimorel (|2000l ) claim that the lack 
of pulse-width evolution in BATSE GRB li g ht cu rves 
favors internal shock models; iKocevski et al.l (|2007f ). in 
the absence of any quantitative consideration of external 
shocks, argue that XRT pulse spreading in Swift GRBs 
favors internal shock models. From our simulations, we 
see that this effect is subtle at best in the unlikely condi- 
tion when the clouds are "uniformly randomly" situated. 
Realistic cloud clustering properties or a narrow 7-ray jet 
could obscure this effect completely. lLazzati fc Pernal 
(|2007f ) derive the time-scale ratio At v /i ~ 0.25 for a 
spherical shell and suggest that the observational data 
indicate a clustering in the value At v /t ~ 0.25 that rules 
out the external shock model, but nowhere consider ob- 
servational biases for the detection of flare profiles with 
larger or smaller time-scale ratios, or whether this time- 
scale ratio could reflect cloud properties. 

As shown in the simulations, STV in the external shock 
model with moderate (~ 10%) efficiency is possible if 
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the blast wave interacts with clouds with thick columns 
and sizes A c ; <C x/Yq, b ecause then the condit ion of "lo- 
cal spherical symmetry" ([Fenimore et al.ll996|) is broken. 
The shortest pulses carrying a significant fraction of the 
total fluence are made by interactions of the blast wave 
with small c louds at angles 9 <C 1/T n to the observer 
line-of-sight (|Dermer fc Mitmanl 120041 ). Efficiency con- 
cerns are reduced if the fraction of energy dissi pated as 
X-ray s and 7 rays is a small fraction of the total (|Dermerl 
f2007f) . 

5. DISCUSSION 

We have critically examined whether external shock 
processes can make short timescale X-ray and 7 ray vari- 
ability in the 7-ray luminous phase and the early af- 
terglow phase, extending to > 10 4 s after a GRB trig- 
ger. The result is simple: external shocks can only make 
the erratic variability observed in GRB light curves if 
the cold fluid shell ejected by the impulsive burst event 
does not spread significantly from its initial character- 
istic width At eX p/c at the time of the GRB explosion, 
where At exp , assumed to be <C 1 s, is the time of the 
energy release in the explosion frame. The question of 
whether lack of shell spreading is a valid assumption, at 
least within the Doppler cone primarily sampled by the 
observer, is an open question but plausibly valid. A cru- 
cial point is that an external shock model for the rapid 
variability is feasible because the GRB blast wave is still 
in its coasting phase, and has not yet reached the adia- 
batic self-similar phase. 

Let us follow the implications if this assumption were 
allowed. 

One way to interpret GRB light curve data is to 
suppose that there is a single impulsive burst of en- 
ergy, and all subsequent phenomenology is a conse- 
quence of the collimated ultra-r elativistic blast wave in- 
teracting with its surrou n dings (jMeszaros fc ReeslUflfll : 
iDermer fc Mitmanl I1999L l2004f ). External shocks are 
the basis for the success that the relativistic blast- 
wave physics model has had in interpreting GRB after - 
glows (jMeszaros fc Reeslfl997t IZhang fc Meszarosl[2003 ) ; 
we also suppose that its success can be translated to 
the prompt phase. The properties of the surrounding 
medium, as implied by extensive evidence on GRB host 
galaxy studies, is determined by the type of GRB pro- 
genitor, which is a massive (> 10 - 20 M Q ) star lacking a 
H envelope whose core collapses to form a Type lb/c SN. 
The SNR and the progenitor stellar wind, including ef - 
fects of a possible binary companion ([Fryer et al.ll2007h . 
form a complicated surrounding s that imprints its struc- 
ture on the GRB lig ht curves (|Dermer fc Mitmanl 1 19991 : 
iWang fc Loel3[2T)Oft) . 

The brief impulsive energy releases in GRB explosions, 
which for neutrino-mediated processes are most efficient 
for ~ ms timescales, are compatible with brief collapse 
events or phase transitions from stellar cores to compact 
objects, for example, from neutron star to black hole 
or neutron star to a q uark star (e.g.. iDrago et al.l 120071 : 
iBerezhiani et al.ll2003f ). This collapse event would follow 
by minutes, hours, days or more a Type lb/c supernova 
that forms a supramassive, rotationally supported neu- 
tron star/magnetar. This picture is a two-step collapse 
mode l, like the supranova model of lVietri fc Steflal (fl9M 
1999), though with a short delay. 



Such a picture, orthogonal in all respects to the col- 
lapsar model except for the primary role of high mass 
stars as long-duration GRB progenitors, has a great deal 
of explanatory and predictive capability. Here we point 
out some of the attractive features of such a picture: 

1. 56 Ni. The direct collapse of the evolved core of a 
massive star to a black hole in the collapsar model, 
avoiding core bounce, requires new routes for the 
production of 56 Ni. This makes the similarity be- 
tween the peak optical rest frame My magnitudes 
(taken as a proxy for 56 Ni production) measured in 
SNe associated with G RBs and Type lb/ c SN e not 
associated with GRBs ( Sod erberg et al.l feOOEl) a re- 
markable coincidence. Two-step collapse processes 
avoid this coincidence by forming 56 Ni through 
neutron-star core bounce. 

2. Energy reservoir. The discov e ry dFrail et al.l 1 2001: 
i Panaitescu fc Kumarl 120011 : iFriedman fc B loom 
2005) of a standard energy reservoir, or at least a 
clustering of beaming-corrected absolute energies 
towards a value = 10 51 ergs, modulo an uncertain 
efficiency to convert explosion energy into 7 rays, 
fixes the energy budget of the GRB explosion. The 
collapse of a rotationally-supported neutron star 
to a black hole, or the transition from a degen- 
erate neutron fluid to quark state, could liberate 
kinetic energies approaching a constant maximum 
energy that may not defy ability or require exces- 
sive assumptions to calculate. If the explosion in- 
volves a neutrino or Poynting-flux mediated inter- 
action, values of collimated electromagnetic energy 
of w 10 51 - 10 52 ergs could be released in a second 
collapse to a black hole. 

3. Complex ci rcumburst medium. A magnetar phase 
(|Usovl 11992ft follows, even if only by minutes, the 
first collapse event, and drives a strong pulsar wind 
to disrupt, thread, and render extremely complex 
the surrounding medium, in particular, the super- 
nova remnant shell produced in the first c ollapse 
event. The work of iKonigl fc Granotl (|2002f ) shows 
the profound impact on the circumburst medium 
that this phase can have for long (month-yr) de- 
lays. 

4. Clean environment. The Usov phase cleans the en- 
vironment except for dense clumps formed at dis- 
tances where the strong radiation pressure can be 
withstood — and these are the clumps that form the 
clouds in the external shock model. The MHD 
winds effectively drive out all surrounding ma- 
terial except for very dense clouds, possibly al- 
lowing for the existence of density contrasts of 
magnitude defined by eq. ([Ml The clouds them- 
selves would be expanding, but on a much longer 
timescale than the period of activity of the system. 
Close to the explosion trigger, the environment of a 
rotationally-supporte d neutron star is baryon-clean 
([Vietri fc Stellalfl998ft . 

5. Explanation for chromati c X-ray and opt i cal af- 
terglow beaming breaks. Panaitcs cu et al.l (|2006l ) 
compare Swift BAT X-ray light curves with opti- 
cal light curves taken from a variety of telescopes 
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for 6 GRBs, finding several instances where the 
X-ray light curves appear to decay chromatically 
with respect to the optical light curves. Possible 
explanations for these observations are evolving e e 
and €b parameters and the addition of the (syn- 
chrotron self-)Compto n for a beamed outflo w (illus- 
trated numerically by iDermer et al.ll2000f) . which 
plausibly explains the light curves of GRB 050922C 
and the lack of X -ray breaks observed with Swift 
(|Sato et all l2007fh The optical light curves of 
GRBs 050607, 050802, and 050713A are too poorly 
sampled to draw definitive conclusions, and the op- 
tical and X-ray light curves in GRB 050319 appear 
to follow each other except for a final high optical 
datum. This leaves only GRB 050401 to explain. 
More and better data is clearly needed. 

6. Explanation for the decay and plateau phases ob- 
served with Swift. Within the context of the exter- 
nal shock model for a GRB taking place within 
a uniform surrounding environment, I recently 
showed that an intense photohadronic phase, effec- 
tively depleting internal GRB blast wave energy, 
could cause t he X-ray light curves to exhibit rapid 
flux declines (|Dermerll2007h . After the discharge, 
the blast wave regains its nonradiative character 
to form the plateaus observed with Swift. Because 
GRBs showing a steep decay phase have signifi- 
cantly decelerated, the blast wave should then have 
entered the adiabatic self-similar phase when only 
broad X-r ay flares can result fr om large density 
contrasts (jNakar fc Granotl[200l . It is interesting 
to note that GRBs showing the most extreme de- 
cay phases, e.g., GRBs 050315, 050416A, 050713B, 
050814, and 050915B, display smooth subsequent 
light curves (|Q'Brien et al.1 [20061 : IChincarini et al.1 
I2007D . though GRB 050916 provides a counterex- 
ample and a challenge to this interpretation. 

7. High energy radiation. Observations of GRBs 
with the GLAST GBM and LAT will monitor the 
Compton components in the spectrum of a GRB. 
Definite correlations between the leptonic syn- 
chrotron and SSC components are expected, which 
behave in stark contrast to photohadronic 7-ray 
components that vary independently of the lower- 
energy lepton synchrotron component. GLAST 
will search for photohadronic emission components 
and, in conjunction with detectors such as IceCube, 
test multiwavelengt h and multi-channel predictions 
(jDermer et al.ll2007f ). 

8. Predictions. With respect to the short-delay supra- 
nova model, the confirming prediction is to glimpse 
the GRB progenitor in its magnetar phase prior to 
the second-step collapse to form the GRB. Analysis 
of radio, optical, X-ray, or 7-ray emissions to ob- 
tain P(t) and P(t) in the brief interval between the 
two collapse events is now only possible with broad 
field-of-view detectors with poor sensitivities. Un- 
fortunately, optical supernova discovery generally 
happens within a few weeks after the supernova 
event, by which time a second collapse would have 
happened. So for the short-delay supranova in- 



tervals consistent with the offset betwe en nearby 
GRB s and their supernova emissions (|Zeh et al.l 
2004), this prediction is not prom ising except to 
test t he standard supranova model (jVietri fc Stellal 
H998h . 

A definitive test of the collapsar model is to find 
a supernova taking place well in advance of the 
GRB, so that the GRB has no detectable super- 
nova emissions. Development of telescopes with 
wide field-of-view optical survey capabilities, like 
Pan-STARRS or SNAP, holds promise to rule out 
the collapsar model for specific GRBs (C. Dermer 
& A. Drago, unpublished, 2006), and could also 
provide a database to search for magnetar activity 
in advance of a GRB. 

6. SUMMARY AND CONCLUSIONS 

A detailed analysis of the interaction between a rela- 
tivistic blast wave and a stationary cloud was performed 
in the limit that that the cloud width A c / <§; x, so that 
the shell density remains effectively constant during the 
interaction. Synchrotron light curves and spectra from 
such interactions were calculated for a range of cloud 
densities, sizes, and locations, for blast-wave coasting 
Lorentz factors To = 100 and 300, and for different blast- 
wave widths. 

In order to produce short timescale variability in the 
external shock model, the shell width parameter 77 must 
be < 1/ro, as demonstrated both analytically and nu- 
merically. If this model assumption is valid, then the ex- 
ternal shock model can explain the generic spectral shape 
of GRB pulses, rapid variability in GRB light curves ob- 
served with BATSE, and the delayed X-ray flares ob- 
served with Swift. The tendency of the GRB pulses or 
X-ray flares to diminish in intensity with time is a conse- 
quence of the expansion of the blast wave. The durations 
of GRB pulses reflect the dimensions of the surround- 
ing thick-columned clouds. The typical distances of the 
clouds that produce 7-ray pulses are w 10 15 - 10 17 cm 
from the GRB explosion center, with cloud sizes ~ 10 12 
- 10 13 cm and densities ~ 10 3 cm -3 . The clouds that 
produce X-ray flares observed with Swift are typically 
found > 10 17 cm from the center of the GRB explosion, 
and have sizes < 10 15 cm and densities ~ 10 3 cm~ 3 . 

This analysis also provides a basis for making accurate 
calculations of light curves and spectra formed by both 
forward and reverse shocks in collisions between relativis- 
tic shells. 
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TABLE 1 
Standard Blast Wave and Cloud 
Parameters 21 





GRB pulse 


X-ray Flare 


Rcdshift z 


1.0 


2.0 


d L (10 28 cm) 


2.02 


4.80 


Eo ergs 


10 53 


10 54 


To 


300 


100 


A 


10 7 cm 




V 


ro\o 




e e 


0.1 




ffi 


io- S io~ 4 




£2 


0.1 




V 


2.5 






0.01 




9i 


0.0 




xo(cm) 


10 16 


10 17 


X2 (cm) 


1.02 x 10 16 


1.02 x 10 17 


n ci (cm" 3 ) 


10 3 





a Forward and reverse shock parameters 
are the same. 



TABLE 2 

Parameters for Simulated Light Curves 3 - 



A B C D E F 



E (ergs) 


10 53 


10 54 






Rcdshift z 


l 








kg 


l 


0.5 






To 


300 






100 


e b 


0.2 




0.002 






0.4 






10 17 


fil(cm) 


10 16 






J?2 (cm) 


10 17 




10 18 




ri(cm) 


3 x 10 12 


10 13 3 x 10 12 


10 13 


3 x 10 13 




0.1 








c 


0.1 









a Only the parameters changed from the previous column are 
listed. b Received photon energy and e p k in units of m e c 2 . 



